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Abstract. In this work we consider weak solutions of the incom- 
pressible 2-D porous media equation. By using the approach of De 
LeUis-Szekelyhidi we prove non-uniqueness for solutions in L°^' in 
CN| space and time. 

O 

Q 1. Introduction 

\Q The incompressible 2-D porous media equation (IPM) is described 

by 

pt + V ■ [vp) = 

where the scalar p{x, t) is the density of the fluid. The incompressible 
velocity fleld 

V ■ = 

is related with the density by the well-known Darcy's law pLj 

p 



< 



K 



-V = -Vp - (0,5(p) 



>■ where p represents the viscosity of the fluid, k is the permeability of 

^ the medium, p is the pressure of the fluid and g is acceleration due to 

fSl gravity. Without lost of generality we will consider p/ k = g = 1. 

^ In this paper we study the weak formulation of this system and we 

CN construct non trivial solutions with p,v E L°°(T^ x [0,T]) from initial 

^ data p(x, 0) = 0. Here is the two dimensional flat torus. 

O We deflne a weak solution of IPM {p,v,p) if Wip,x, ^i, ^2 G 

> C^{[0, T) X T^) with A = (Ai, A2) the following identities hold 

(1) J J p{x,t) {dt(f{x,t) + v{x,t) ■Vip{x,t)) dxdt 

T2 

+ J poix)Lp{x,0)dx = 0, 

J2 



(2) yy„(M).vx(MM.*^o. 

T2 
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1 

(3) j j {v{x,t) + Vp{x,t) + {0,p{x,t))) ■ \{x,t)dxdt = 0. 

T2 

For initial data in the Sobolev class if*(T^) {s > 2) there is local- 
existence and uniqueness of solutions in a classical sense and global 
existence is an open problem [1]. It is known the existence of weak 
solutions, where the motion takes place in the interface between flu- 
ids with different constant densities, modeling the contour dynamics 
Muskat problem [3] . The existence of weak solutions for general initial 
data is not known. In this context we emphasize that the solutions we 
construct satisfy 

limsup = +00 



for any s > (see Remark 5.1). 



From Darcy's law and the incompressibility of the fluid we can write 
the velocity as a singular integral operators with respect to the density 
as follows 

v{x,t) = PV [ n{x-y)p{y,t)dy -\{{],p{x)), x G M^ 



where the kernel is a Calderon-Zygmund type 



27r V kl^ 

The integral operator is defined in the Fourier side by 

^(0 = (^,-^)p(0- 

This system is analogous to the 2-D surface Quasi-geostrophic equa- 
tion (SQG) [Ij, in the sense that is an active scalar that evolves by 
a nonlocal incompressible velocity given by singular integral opera- 
tors. It follows that, for Besov spaces, if the weak solution p is in 
L3([0,T] X ^3'°°) with s> I then the norm of p is conserved |27j . 
This result frames IPM in the theory of Onsager's conjecture for weak 
solutions of 3-D Euler equations [2], [16]. However there is an extra 
cancelation, for SQG, due to the symmetry of the velocity given by 

that provides global existence for weak solution with initial data in 
L^(T^) [20] • Furthermore, one can find a substantial difference be- 
tween both systems for weak solutions of constant p in complementary 
domains, denoted in the literature as patches [I3]. For IPM the Muskat 
problem presents instabilities of Kelvin-Helmholtz's type [3J and there 
is no instabilities for SQG ([2T].[8]). 
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The first results of non-uniqueness for the incompressible Euler equa- 
tions are due to Scheffer [22] and Shnirelman [23] where the velocity 
field is compacted supported in space and time with infinite energy. 
The method of the proof we use in this paper is based on understand- 
ing the equation as a differential inclusion in the spirit of Tartar [ZS1I2S]- 
In a ground breaking recent paper De Lellis-Szekelyhidi [6j showed that 
with this point of view the modern methods for solving differential in- 
clusions could be reinterpreted and adapted to construct wild solutions 
to the Euler equation with finite energy. Our plan was to investigate 
the scope of this approach in the context of the porous media equation. 
However it turned out that the situation is different and we have to 
take different routes at several places which might be of interest for 
the theory of differential inclusions. We describe them shortly. Firstly, 
using the terminology of this area that will be recovered in the first 
section, the special role of the direction of gravity yields certain lack of 
symmetry in the wave cone A. Moreover the set K describing the non 
linear constraint belongs to the zero set of a A convex function. Then 
it follows that, opposite to Euler, the A convex hull does not agree with 
the convex hull and more relevant K C dK^. This is an obstruction 
for the available versions of convex integration, the ones based on Baire 
category ^ [TUl US] and direct constructions [HI HB E] • Thus we do 
not see anyway to use the method to produce solutions to the equation 
with constant pressure and velocity for some time. Surprisingly, an 
easy argument shows that the states in int(K'^) can still be used to 
produce periodic weak solutions starting with p = f = 0. This leaves 
the difficulty of choosing a proper subset K G K with sufficient large 
A hull. In general the computation of A hull might be rather compli- 
cated. In this work we argue differently to suggest a more systematic 
method: Instead of fixing a set and computing the hull, we pick a rea- 
sonable matrix A and compute {A + A) (1 K. Then by [10, Corollary 
4.19] it is enough to find a set ^ C {A + K) r\ K such that A e 
to find what are called degenerate T4 configurations supported in K'^. 
To our knowledge this is the first time they are used to produce exact 
solutions. Then, we are able to choose the set K carefully so that the 
construction is stable and allow us to solve the inclusion. Our proof of 
this later fact is related, to the known ways of solving differential inclu- 
sions based on T4 configuration [2H |lTl |10l HI] but we have made an 
effort to keep the whole paper self contained. In particular we do not 
appeal to the general theory of laminates [TSl HHl [E] to minimize the 
amount of new terminology and make the construction more explicit. 

The paper is organized as follows. In Section 2 we review the De 
Lellis-Szekelyhidi approach to obtain weak solutions to a nonlinear 
PDE . Section 3 is devoted to put IPM in this framework, the compu- 
tation of the wave cone and the corresponding potentials. Section 4 is 
devoted to the geometry part of the construction and Section 5 to the 
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analytic one. In this last section we include some remarks about the 
nature of our solution as well as why this approach seems to fail for 
constructing weak solutions compactly supported in space and time. 

2. General strategy and notation 

The strategy to construct wild solutions to a non linear PDE , under- 
standing it as a differential inclusion, builds upon the following main 
steps. Consider the equation^ 

P{u) = 0, u:W ^ W^, 

where P is a general differential operator. 

Tartar framework: Decompose the non linear PDE into a lin- 
ear equation (Conservation law) C{u) =0 and a pointwise constraint 
u{x) G K. Possibly this requires to introduce new state variables u to 
linearize the non linearities 



P(u) = 



C{u) = 

u{x) e K 

Specially easy to bring to this framework are equations of the form 



C{Fi{u)) = 

where Fj : M™ M i = 1,...,/ are some functions possibly non 
linear. Then one considers new states variables Ui = Fi{u), which 
clearly satisfy the conservation law. Then the pointwise constraint 
amounts to 



Qi = Fiiu). 

Wave Cone: As discovered by Tartar [25] in his development of 
compensated compactness to such conservation law one can associate 
a corresponding wave cone A^, defined as follows: 

A^ = {A e X M' : 3^ such that for z{x) = Ah{x ■ , C{z) = 0}, 

where /i : M — M is an arbitrary function. That is, the wave cone 
corresponds to one dimensional solution to the conservation law. 
Potential: li A e Ac we can build ^ : ^ M"' x M' such that 

C{z) = and C 

where = {tA : t G M} is the line passing through A. The wild 
solutions are made by adding one dimensional oscillating functions in 
different directions A. For that it is needed to localize the waves. This 
is done by finding a suitable potential, i.e. a differential operator D 
such that 

C{D) = 0, 

^Here n and m are suitable spaces dimensions that we keep general at this point. 
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and the plane wave solutions to the conservation law can be obtained 
from the potential. 

Pointwise constraint: In the pointwise constraint one adds the 
additional features desired for the weak solution (fixed energy, fixed 
values, etc). This describes a new set K G K. Then the solution is 
obtained by two additional steps. 

• Find a bigger set U so that there exists at least a solution to 
the relaxed system 



with the proper boundary conditions. In practice the set U C 
mt{K^), the so-called A convex hull of K. 
• Then one fixes a domain Q compact in space and time and 
in properly chosen subdomains of Q one adds infinitely many 
small perturbations to find a function Uoo € K a.e (x, t) & Q 
and Uoo = Uq outside ofQ. If Uq is a solution to the conservation 
law, so is Uoo- Notice that then Uq & K and to be sure that 
Uqo 7^ Uq, we need that Uq ^ K. 

There is no standard way to produce the weak solutions from 
and it is not even always possible. In (for Euler equations) it is 
enough that (0,0) eU = K'^, which is open in the spaces where the 
images of the potential lie and that for every state U eU, Au E A such 
that the segment which end points u± A eU and 

\A\ > Cdist{u,K). 

As discussed in the introduction, this strategy fails in our context 
and needs to be modified. 



Our task below is to adapt section 2 to the IPM system. 

3.1. Porous media equation in the Tartar framework. Given 
functions (p, v, q) we define the differential operator C by 



UqEU 



3. Porous media equation 



Here we start by recalling the 2-D IPM system 

[dtp + V ■ ivp) = 





'dtp + V ■ g 



V ■ V 



(5) 



C{p,q,v) = lCm\{v + {Q,p)) 
^/t2 P 
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Proposition 3.1. A pair {p,v) G L^{T^ x [0,T],M x M^) ^5 ^ ^^ak 
solution to IPMifand only if we find u ^ {p,v,q) e L°°{T'^ x[0,T],Rx 
X R2) such that, 

C{u) = 
pv 

Proof. The two first equations are obvious. The third follows from 
the Hodge decomposition in T^. Namely consider a vector field / e 
L^(T^) with fdx = and curl / = 0. Wc claim that the Hodge 

decomposition implies that / is a gradient field. Namely, every field 
/ e L^(T^) on the torus has a unique orthogonal decomposition 

+ " + 

such that divw = 0. But curl/ = = curlw. Now we have that 
divw = and curlw = 0, then for example by using the formula 
curl curl gf = — Agi + Vdivgi one finds that Aw = and which yields 
w = 0. Thus, 

□ 

It is convenient to write the differential part of £ in a matrix form. 
For that to each state u = {p,v,q) we associate a matrix value function 

U{u) : R X R2 X R2 ^ M=^^^ 

by 

(~V2-p Vi, 0\ 
Qi Q2 p) 

and the subspace U C M^^^ ^y U = ^/(R^ x R^ x R). 
The following lemma is straightforward. 

Lemma 3.1. Let U = U{p,v,q) G L^(T^ x R) with pdx = 

fj2 vdx = 0. Then the following statement holds in the sense of 
distributions: 

C{u) = Div{U) = 0, 

where Div{aij) — djUij 
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3.2. The wave cone. We denote by A the wave cone related to (IPM). 

The reason to write the states in matrix form is that the wave cone 
is particularly easy to characterize. Namely, we want a plane wave 
z{x) — Ah{x • ^) to satisfy 

Div{z) = 0. 

We obtain that for each i, h'{x ■ ^)Qaij = 0, that is 

A{0 = 0. 

In other words A G A if and only if ^4 is a singular matrix. Thus 
yl e A n U if 

—p\V2 + Vi + PV2) = <^=^ p — or p\v\ — 0. 

V2 



Observe that fixed p this is 

vl + vj + pV2 =vf + {V2 










HI 












0/ 



2^ 4 

Thus fixed p the corresponding i'eS'((0,— ^),^). 

3.3. Potentials. The aim of this section is to find a differential oper- 
ator D such that C{D) = 0. We define for ip,(p : x R ^ R, 

dxiX2^ -dx^xii^ 

-dtxiip - dx2<P - dtx2 ip + dx^(p AV'^ 
Or more compactly written 

-dx2X2^ dxiX2fP 

dxiX2fp -dxix^ip 

-dtxii^-dx2V -dtx2'ip + dxi^p Aip^ 

Clearly D{ip, -0) G U. 

Lemma 3.2. Let Q G x R be a bounded open set, (p G W^'°°{9,) 
and V' e W^'°°{fl). Then it holds that 

jCD{ip, -0) = m Q 

in the sense of distributions. 

Proof. In our notation we need to obtain that 

Div{D{(f,ilj)) = 

distributionally. The claim follows since distributional partial deriva- 
tives commute. Namely, it holds that 

I* 9x2X2Xl'4^ ~l~ 9x^X2X2^ — ^1 

II: dxiX2Xi'4^ 9xixiX2^ — 0) 
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III: {-dtxix^i^ - dtx2X2i^ + dt/^i)) + {-dx2xi<^ + d^^^i^) = 0, 

as desired. Finally we observe that since the flat torus has no bound- 
ary then 



Remark 3.1. Observe that we are solving Atp = p and then using the 
PDE to define the rest. 

3.4. Plane waves with potentials. In this section we show that for 
any direction A with p 7^ we can obtain a suitable potential. It will 
be more convenient for us to work with saw-tooth functions instead 
of trigonometric functions. The following proposition is related to [101 
Proposition 3.4]. 

Lemma 3.3. Let U e A with pv ^ 0, < X < l,e > 0,n C x R. 

Then there exists a sequence ujy : fl —>■ V of piecewise smooth functions 
such that 




and thus CD{(p,ip) = 0. 



□ 



• C{un) = 0, 

• sup dist{uj\f{x,t),[—{l — \)U,\U])<e, 



• \(x,t) e n : UNix,t) = (-1 - \)U\ > A(l - e), 

• \{x,t) e Q : UN{x,t) = XU\ > (1 - A)(l - e), 

• UN ^0 %n L°°(T2 X M). 



Proof Let S G VT^.ooj^^^^) ^^^^ 




and 




Then 



DidipN,iJN) = S"iNi 




Xi + 




V2\x2+Ct)) 



/ -V2 - P 



Vl 

V2 








Now the matrix 
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has determinant —{vf + V2)/\v2\ ^ and hence defines a bijection of 
M^. Thus for any qi, q2 we can choose (c, d) so that 

- dVWli -ca/H - d '"J— ) = (gi, q2)- 

Therefore we have that for any direction in A with pv ^ we can 
obtain a suitable potential. 

Next we choose an appropriate function S. We consider the Lipschitz 
1-periodic functions S", s : M — M such that S{0) = 0, S' = s, s'{x) = 
1 - A - X[A/2,i-A/2](a;) for x G [0, 1]. We have that s(l/2) = 0, s'(l/2 + 
x) = s'(l/2 - x) and s(l/2 + x) = s(l/2 - x) if \x\ < 1/2. 




We then localize the wave using the potential in the standard way 
i.e. we define a test functions Ce' such that 

|Ce'|<l, Ce' = 1 on 5w(0), supp(C') C fil(O). 

Let Q = -Bi(O) and for a suitable e' we define = D{Q{ipiy,ipN)). 
For a general Q consider disjoint balls Br^.{xk,tk) such that 



\n \ UkBr,{Xk,tk)\ > 1 



and finally we define 



UN{X,t)=U^[ , ) on X e Bri^{Xk,tk). 



□ 
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4. The construction 

4.1. Geometric setup. In this section we identify U with M x x M^. 
For A eV we will use coordinates (p, w, z) with p e M and w,z E M^. 
To manipulate the A cone it will be helpful to introduce the following 
notation. Along this section we denote by 



Sp = 5((0, 



2 2'' 
Similarly, 

The set K = {p,w,pw) C U defines our pointwise constraint. The 
strategy in Section 2 requires to find Xq e {K\K) nint(i^'^). However 



Proposition 5.1 shows that this is not possible. Instead we will find 
states of the form (0, 0, z) G int(i^'^). In the next section it will be 
shown that this is enough to produce a weak solution. 

We choose a suitable K C K. For that we need to introduce the T4 
configuration. 

Definition 4.1. We say that C G U is the center of T4 configuration 
(degenerate) if there exists {Ti{C)}f^i G U such that 

a) c G mmur, 

b) C-Ti{C) G A. 

We say that A belongs to the T4 configuration if A G [C,Ti{C)] for 
some i. In this case we also set Ti{A) = Ti{C). 

T4 has been the key pieces of subtle versions of convex integration 
112]. In this degenerate form they appear for the first time in |lUj . 

Lemma 4.1. Let zq ^ (0) "2^) ^ -^i- Then there exists S = 
S{zo) such that every A G U with \A — {0,0,zq)\ < 6 is the 
center of mass of a T4 configuration in K. The Ti{A) are of 
the form (1, a;i, si), (1, 0:2, 0:2), (-1, yi, -yi), (-1, ?/2, -^2) G K, where 
{xi^X2iyiiy2){A) : ^ are differentiahle submersions. Moreover 
A G int{Ti{A)Y. 

Proof. First we observe that by the definition of the cone A 
(p, w, z) — (1, X, x) G A '^=> X — w E Si^p 

and 

(p, w, z) - (-1, y, -y) G A y-w e 

if in addition we have 

ip,w,z) G int({(l,xi,xi), (1,X2,X2), (-1,1/1, -yi), (-1,1/2,-1/2)}'') 
the lemma is proved. 
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First, we notice that for each (p, z) with |p| < 1 there exists 
different iable functions x, ?/ : M x x ^ such that 

Namely for t = and 

/ N Z + W W - Z 

x{A) = x(p, w, z) = — — , y{A) = y{p, w, z) = , 

1 +p 1 -p 

it holds that 

ip,w, z) = (t(l) + (1 - t){-l),tx + (1 - t)iy),tx + (1 - t){-y)). 

Now we consider the balls = w + -Bi-p, By = w + -B-i-p with 
centers ay. 

Thus, if 5 < 5{zq) is sufficiently small it follows that 

(6) X = e Bx\ax, y = ^ _ ^ ^ By\ a^ 

In this case we can use the convexity of Euclidean balls to show that 
x{A) G [xi, X2], xi, X2 E w + Si-p and similarly for y{A). Namely if 
\x — a\ < r then 



X = -(1 - -) (a - r- -) +-(1 + -) (a + r- 



2 r \x — a\ 2 r \x — a\ 



X2 
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The same argument provides us the mapping yi{A),y2{A) Thus we 
declare 

Ti{A) = (l,Xi,Xi),T2(A) = (1,X2,X2), 
T,{A) = (-l,y,,-y,),T,{A) = (-1, 1/2,^2). 

We have shown the existence of {Xi}f^i such that < Aj < 1, 
Ylt=i -^i = 1 and 



1=1 



Since the Ti{A) are hnearly independent and Aj 7^ for i = 1, 2, 3,4 
we arrive to the desired claim 

A G int(U,T,(A))^ 

We turn to the regularity of the mappings Ti{A) : W ^ M^. It 
is enough to consider the case Ti. Notice that from the convexity 
argument follows 

(7) = a{w,p) = w + {1 - p){0,^) and r = 1 - p. 

Thus, both are differentiable functions of A. Since x ^ the same is 
true for xi,X2,yi and 2/2- It remains to see that they are submersions. 
It is enough to argue for the function xi. We write xi = Fi o F2 with 
Fi : M — >• M^, Fi{r,z,w) = {rx,ax,x) which are defined in (7][6) and 
F2 : — >• is given by F2(r, a,x) = a — r-^^^. For small 5 Fi is a 
submersion and F2 is always a submersion. □ 

Remark 4.1. It is instructive to realize that by definition if {p,w,z) 
belongs to n A then Z2 < and thus (0, 0, 0) ^ int(i^' n A)^ This 
shows that the above strategy never could work with z = (compare 



with Proposition 5.1) 



We will choose a set U adapted to the T4 configuration. 

Definition 4.2. Let K cV. The first Lambda convex hull of K, K^'^ 
is defined by 

K^'^ = {AeV : A = tB+{l-t)C, B,C e K, B-C e A,t e [0,1]}. 
Definition 4.3. Let z 7^ (0, ^) G Bi and the corresponding S{z) from 



lemma 4.1 , We denote by G U the euclidean ball centered in (0, 0, z) 



and with radius S{z). Furthermore we declare 

= uti{B,UT,(B,)}i'^ and = uti{r,(B,)} G iT. 
Lemma 4.2. The set Uz\K is open in U. 



POROUS MEDIA EQUATION 



13 



Proof. Let C e U with |C| < e and At eUz\ K. Our task is to find 
an e so that Ct = At + C eU^\K . By definition At = tA + {l- t)XA, 
where A e B^, < t < 1 and G Ti(B^) for some i = 1,2,3,4. 
Moreover it holds that A — € A. 

It will be enough to show that if e is small enough, then there exists 
X e Ti(M,) such that 



Ct-XeA and X + ^{Q - X) e M,. 

We use coordinates A = {p,w,z), Xa = (l,^^,^^). At = {pt,Wt,Zt), 
C = {pc, wc, zc)- Since Xa — A e A it holds that Xa — At & K as well. 
Thus, by the definition of A, xa — Wt & Si-p^. Therefore there exists ^ 
with 1^1 = 1 such that 

XA = wt + {i-pt)m~^)+^]. 

Next recall that Ct — G A if x — — wc G 5'i_(p^_|_p^). Set 

X = Wt + Wc + (1 — — pc')((0, ^) + notice that 

\x - xa\ < \wc\ + 2\pc\ < 3\C\ < 3e. 



The function x : A — > x{A) is a submersion by Lemma 4.1 and 
thus ^(B^) is an open subset of M^. Hence, by choosing 6e < 
dist(x^, ^^(B^)), we obtain that (l,x,x) G Ti(Mz). 

Now since A = Xa + ~ ^a) "we have 

\X+-^iCt-X)-A\ = \X~T,{A)\ + j\Ct-At\ + \X-T.,iA)\ 

< (3+^).. 

Then if e < |t ■dist(A, 9B^) it follows that Ac = X +\{Ct-X) eM^ 
and 

Ct = tAc + {l-t)X eu,. 

A final choice 8e < min{t ■ dist(y4, SB^), dist(xA, ^^(B^))} yields the 
claim. □ 



5. From geometric structure to weak solutions 

There are many ways to pass from solutions in U to solutions in K 
depending on the geometry of the sets. We have followed to present an 
approach which in some sense axiomatize the arguments in [B]. However 
it is by no means the only argument and for example one can verify 
that our sets K and U verify the conditions in [lOi Proposition 4.42]. 
The argument is easier to explain via the following definition which is 
related to the notion of stability near K given in [lOl Definition 3.15] 
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Definition 5.1 (Analytic Perturbation Property, APP). Tlie sets K C 
U have the analytic perturbation property if for each A & U and for 
every domain i7 C x M there exists a sequence Zj of piecewise 
smooth function such that 

i) = 0. 

ii) Zj is supported in VL. 

iii) A + Zj & U a.e. 

iv) Zj 0. 

v) There exists c > such that J^\Zj\ > c ■ distK{A)\Q\. 

Next we prove that our sets enjoy the analytic perturbation property. 
The reader that is familiar with convex integration will realize that in 
practice our arguments are related to the concept of in approximation 
jl4t [9] for A convexity (based on degenerated T4 configurations). 

Lemma 5.1. The sets andUz have the APP property. 

Proof. We will obtain v) as a consequence of 
v) There exists ci, cq > such that 

|{(x,t) G Vt : \Zj{x,t)\ > cidist^(A)}| > cq\Vl\. 

Let C G B^. By lemma |4.1| is the center of a T4 configuration 
{Ti{C)}1^i. Now observe that for s small enough C belongs also 
to a T4 configuration supported in Tj ,,(C) = sC + {1 — s)Ti{C) G 
'vaiKz^s- Namely, since Ti^siC) Ti{C) as s — * it follows that 
C G int({Tj(C)})'^ implies that C G {Tj s(C)}'^. The other properties 
are obvious. 

Let Ti^s{C) be as above. First notice that since \C — Ti{C)\ > 
dist{C,Kz) by continuity we can choose sq such that for s > sq it 
holds that 

(8) 2|C-T,,,(C)| >dist(Ai^.). 

Since C G {Ti^s{C)Y, there exist Aj > with ^ Aj = 1 such that 
C = J2i ^i'^i,s{C). Since the T4 configuration is degenerate, if we set 
Cl = C + e Ei=i(^i,s(C) - C) and Tl{C) = T,^s{C) + Q it follows that 

with Tfj^i{C) — C- G A, i.e T/ is a non degenerate T4 (We drop the 
dependence of s for simplicity). 

Fix a sequence r]j such that n^j^(l — rjr) > ^ and e = > 0. We 
will obtain the sequence {Zj}'^^ recursively. We claim that given Zj 
there exists ^j+i such that 

i) \{{x,t) : |Z,+i(x)| > |dist^(C)|}| > \{{x,t) : \Z,{x)\ > 
|dist^(C)|}| + (l-e)%i. 
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ii) There exists i = i{j) such that 

n^,=,(l - vr){l - ey < \{{x,t) : C + Z,^,{x) = Cn\. 

iii) Zj 0. 

iv) C + Zj eU,. 

The proof follows by induction with Zq = 0. We start by considering 
the open set Qj = {{x,t) e Q : C + Zj = Q}. Since Q = (1 - 



i^)Ci_i + j^Tl{C) we apply the Lemma 3.3 in ilj with r] < r/j+i < r]o 
to obtain a new sequence {Zj^k}- 

. \{{x, t) e : + Q = TaC)}\ > (^)(1 - 

. t) e : Z,,^ + Q = Q_ J| > (^)(1 - r^,+i)|fi,|. 

. sup dist{Z,,,{x,t) + Cl[T,^{C),CU])<{l-7^,). 

Now set Zj^i^k = + and e = Let us verify properties 
i) and ii) . Notice that if {x,t) e flj and ^j,i!c(a;) + = Tl{C), then 
= T^{C) — C and thus by (fsl), for small e, we get 



For i) we have 



> ^distA-(C). 



|{(x,t) : > distK(C)}| > \{ix,t) G : > dist;,(C)}| 

+ \{ix,t)enr-z,A^,t) + c^ = T:{c)}\ 

> \{{x,t) en\nj: \Zj{x,t)\ > dist;^(C)}| 
+ e{l - 7]j+i)\nj\ 

> \{{x,t) en\nj: \Zj{x,t)\ > distK{C)}\ 
+ eUll\il-rir)il-ey\n\ 

>ls{i-sy\n\. 

As for ii), we have 
\{{x,t)eQ:Z^+,Ax,t) + C = CU}\ 

> \{{x,t)eQ, : Z,{x,t) + Z^Ax,t) + C = Cl,\ > {l-e){l-r],+,)\n,\ 

>{l-e)^^'Ull\{l-r^,^m- 

Next, property iv) follows because the segments [T/(C),Q_;^] are 
compact subsets of Uz and by a diagonal argument we choose a subse- 
quence such that 
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It just remains to show that Zj satisfies property v in the definition 
of analytic perturbation property (APP). Indeed by property i) in the 
definition of Z^ yields 



G ^ : > ^dist;,(C)}| = ^(1 - tyt\ = - 

k=l 

> -\Q\ 

for j large enough. Thus for C G (APP) holds. 

Let A eU,\ M,. By definition A = sC + {1 - s)Xc with C G B^, 
Xc G Ti(Mz), < s < 1 and C — Xc G A. Consider the sequence 



{Zj} given by Lemma 3.3 If s < | {Zj} is the sequence required by 



the (APP) property. Hence we may assume that s > |, obtaining that 
\{x eQ: A + Zj{x) = C}\ > |. Let Qj = {{x,t) G Q : A + Zj{x,t) = 
C}. Since C G B^ we can use the argument above for the domain Qj. 
We are provided with a sequence {2'^} such that Z^ ^ 0, Zi^ + C 
and 

\{{x,t) G n,- : \Zu{x,t)\ > ^distK(C)}| > > 

A proper subsequence Wj = Zj + Zk{j) satisfies all the desired prop- 
erties in the definition of (APP). As a matter of fact just property v 
needs verification. Notice that for C G B^ and X G Tj(B^) it holds that 

l-6<\C-X\<l + 6 + max \Ti{C)\ = M{6). 
Thus distx(C') > 1 — 6. Next we notice that 

(9) disti^(A) <\A- Xc\ < s\C - Xc\ < M{6). 

Let {x,t) G Qj be such that \Zk{x,t)\ > ^distK{C). Then 



\W,{x,t)\>^distKiC)-\C-A\ 



(10) >^{l-6)-il-s)\C-Xc\ 



Hence if we declare ci{5) = g^^, and we put together (9) and (10) 



we obtain 



Ci{5)distKiA) < \Wj{x,t)\ 
for every {x,t) G Qj with \Zk{x,t)\ > ^distxiC). We arrive to the 



estimate 
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\{{x,t) e n : \Wj{x,t)\ > ci{S)distK{A)\} > -\n\. 
Finally (v') holds with Ci(5) = Ci and Cq = ^. 

□ 

Next we show that if a set has a perturbation property there are 
many solutions for our inclusion. We formulate the existence theorem 
in the following way: 

Theorem 5.1. Let K,U be two sets that satisfy the Analytic Pertur- 
bation Property and C x M a bounded open set and A &U. Then 
there are infinitely many solutions U such that 

(1) £{U) = 

(2) U{x,t) e K a.e {x,t) e O 

(3) ulx, t)^A for {x,t) eT'^ xR\n 

We shall prove the theorem using Baire category since it yields in- 
finitely many solutions. In addition we will provide a direct argument 
to construct such solutions since it helps to understand its nature. 

Definition 5.2 (The space of subsolutions Xq). Let A e U. We say 
that U e L~(T2 X M) belongs to Xq if 

(I) U is piecewise smooth. (Regularity) 

(II) For t < and t > T U{x,t) = A. (Boundary conditions) 

(III) £{U) = 0. (Conservation law) 

(IV) U{x, t) eU a.e {x, i) e x R. (Relaxed inclusion) 

We endow Xq with the L°° weak star topology and declare X — Xq. 
If a set has the APP property then it holds the perturbation lemma. 

Lemma 5.2 (Perturbation Lemma). There exists C > such that for 
every Uq G Xq we can find a sequence Uk G Xq such that Uk ^ U but 

f \Uk - Uof >C [ disf^iUo). 
Jt^xr Jt^x[o,t] 

Proof. The sequence Uk will differ from Uq only on the compact set 
X [0,T]. Since Uo is picccwisc smooth we can find a finite family of 
pairwise disjoint balls Bj x [0, T] such that 

(11) / distK{Uo)<2y2distK{Uo{xj,tj)\Bj\ 

Vt2x[0,T] 

Next we use that UQ{xj,tj) G to apply the APP property in the 
domain Bj. Let us called the corresponding perturbation sequence by 
{Zj,k}'kLi supported in Bj. Thus we consider the sequence 

Uk — Uq-]- ^ Zj^k, 

which satisfies that Uk &U (a small continuity argument is needed). 
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Then 



/ \Uk-Uo\=y2 \Zj^k\>cy2^^^^K{Uo{xj,tj))\Bj 

J ■' 

distil (t/o) 



where the last inequahty follows from (11). 



□ 

Next we recall how to obtain a solution U{x, t) E K a.e. from the 
perturbation lemma and the non emptiness of Xq. 

The solution U will be the strong limit of a sequence Uk which is 
obtained from the perturbation lemma. Hence strong convergence 
implies that U is in K. The Baire category argument yields the 
existence of such U as point of continuity of the Identity map and 
then a sequence Uk which converge to it. The direct argument builds 
the sequence Uk iteratively and then show that indeed it converges 
strongly to some U . 

Baire Category: 

We consider the Identity as a Baire-one map from (X, weak*) to 
(X,LHT2 X M)) (see P Lemma 4.5]). 

Lemma 5.3. Let U E X he a point of continuity of the identity. Then 
U eK a.e {x,t) eQx [0,T]. 

Suppose that U G X is a point of continuity of the identity. By 
definition of X, there exists {Uj} G Xq such that 

Uj-U 

in the weak star topology. As m is a point of continuity of the identity 
it holds that 



lim [ \Uj -U\=0 



and thus, 

lim / distKiUj)dx = / distx(f/). 

j^oo Jj2 X [o,T] Jj^ X [o,r] 

For each Uj we use a perturbation lemma. The corresponding se- 
quence is indexed by Uj^k- By a diagonal argument we can choose a 
sequence Uj^k(j) such that 

and since U is a point of continuity yields 



lim / |f/.^.(,) - [/| = 



Now 
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/ 



distA'(f/) = lim / distii-(f/j) < C lim / \Uj — Uj^k{j)\ 

T2x[0,T] Jj^x[0,T] ^^°°Jj^xRt 

P. Lemma 

<C\mv I {\Uj-U\ + \U -Ujk{j)\) = ^- 

Proof of Theorem 5.1: It is a well known fact that the points of 
continuity of a Baire one map are a set of second category (see [T7j. 
Therefore there are infinitely many U & X such that U & K a.e. and 
{x,t) G Q. This is equivalent to property (II) in Theorem 5.1. Since 
(I) and (III) hold for any element in Xq they hold as well for any U & X. 

Direct Construction: 

Lemma 5.4. To each function Uq G Xq we can associate a sequence 
{?7j}°^o ^ ^0 such that there exists Uoo G X satisfying 

• U^eK a.e. {x,t) G x [0,T] 

• \imj^^J.^,^^\Uj-U^\ =0 

Proof. We construct directly a function Uoo in X and a sequence Uk G 
X which converges strongly to U as follows. We start with Uq in the 
lemma and obtain with f/^+i from Uk by the perturbation lemma. To 
each Uk we associate a moUifier rjj^k) such that 



(12) / \Uk - vm * Uk\ < 

Then we apply the perturbation lemma to the function U to obtain 
a corresponding sequence Us E X. By the weak convergence property 
we can choose s{k) so that for each i < k 



(13) 



T2xIR 



\Uk - Us{k) * Vj{i))\ < 2" 



and finally we declare f/^+i = Us(k)- Properties (12) and (13) yield 
strong convergence. Now 



/ distx(f/) = lim / dist 



K{Uk 



< C lim 



p. Lemma 



\Uk-U 



k+l\ 



T2xn 



□ 



Theorem 5.2. For every T > there exists infinitely many non trivial 
weak solutions {p,v) G L°°{T'^ x [0,T]) to the 2D IP M system (Q such 
that p{x, 0) = 0. 



20 



DIEGO CORDOBA, DANIEL FARACO AND FRANCISCO GANCEDO 



Proof. We fix 2 7^ (0,—^) G Bi and the corresponding B^, Uz and 



Kz given by definition 4^ By Lemma 5A_ Uz and Kz enjoy the APP 
property. Thus we can apply Theorem 5.1 to these sets, the domain Vt 
equals x (0,T) and A equals (0,0, z). 

We claim that each of the functions U = (p, v,q) E X with U (x, t) G 
Kz a.e {x,t) G X [0,T] from Theorem 5.1. yields a weak solution 
(p, v) to the IPM system with initial values p{x, 0) = 0. Indeed from 
the fact that £(p, v,q) = it follows that div(t>) = and that Curl(t> + 
(0,p)) = 0. Thus there exists a function p such that Vp = v + (0,p). 
The only delicate issue corresponds to the equation of conservation of 
mass where we see the nonlinearity. By definition we have that for 
ip G C^[T^ X M] 

= / (dttpp + V^P-q). 

Notice that G C^[T^ x [0, T)] can be extended to G C^[T^ x M]. 
Next we observe that q{x,t) = z for a.e. {x,t) G x M \ [0,T). Thus 
for a.e. t G M \ [0, T) it holds that 

(14) / Vijj{x,t) ■ q{x,t)dx = / Vip{x,t) ■ zdx = 0. 

Now we plug in (|5]). Then (14) together with the fact p(x, t) = for 
a.e. {x, t) G T2 X M \ [0, T) yields 

T T 

= j j{dtW + Vv-q) = j j{dtVp + Vv-pv), 

T2 T2 

where in the last equality we have used that for (x, t) G x [0, T], g = 
pv. Hence (p, v) is a weak solution to the IPM system according to ([T]), 
(§, (§. □ 

Remark 5.1. The boundary values (or initial data) are attained in 
the weak sense. Namely we use the fact (see [7]) that a solution to 
^{U) = can be redefined in a set of times of measure zero in a way 
that the map 

t ^ v{x, t)ip 



is continuous. Since for negative t J^j v{x, t)(p = we obtain the desired 

lim / v(x, t)ip = 

and we can argue similarly for the density. From this point of 
view it is not surprising that lim^^o p(2;) ^) = = lim^^o ^(^^^ ^) but 
limt_^o p{x, t)v{x, t) = z since they are weak limits and weak limits 
do not commute with products. In fact from this we deduce that our 
solutions do not attain the boundary values in the strong sense and 
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hence by the Frechet Kolmogorov theorem they are highly irregular. 
Another additional feature of the weak solution is that t)\ = 1 for 
a.e. {x,t) G X (0,T) which prevents the pointwise convergence to 
the initial data. 

We conclude with a proposition pointing out to the fact that this 
method does not seem to yield a solution compact in space and time. 
For that we need some more standard terminology. We say that a 
function / : U ^ M is A convex if for A,B,C eV such that A-B e A 
the function g{t) : f{C + t{A — B)) is convex. Let K dJ]. Then the 
Lambda convex hull of K , 

K^ = {Ae\]: f{A) < sup /(C), C eKJ A convex }. 

In all the methods to solve differential inclusions the canonical choice 
for the set U is the interior of K^. However in our case the set K can 
not be there. 

Proposition 5.1. The set K C dK^ 

Proof. The proof relies on the fact that the function f{A) = Z2 — pv2 
is A convex. In fact it is easy to see that detA is A linear, thus — detA 
is A convex. Hence 

If P — det = —pv2 

A 

is A convex and / is A convex. Moreover f{A) = for all A & K and 
hence C f-\-oo, 0] and K C df-\-oo, 0] n iT^ C dK^. □ 

Remark 5.2. The above argument also shows that if G (or 
any other element in K) the elements of the splitting sequence in the 
laminate [TU] should belong to /~^(0). 
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